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We propose new types of density dependent contact pairing interaction which reproduce the 
pairing gaps in symmetric and neutron matter obtained by a microscopic treatment based on the 
nucleon-nucleon interaction. These interactions are able to simulate the pairing gaps of either the 
bare interaction or the interaction screened by the medium polarization effects. It is shown that the 
medium polarization effects cannot be cast into the density power law function usually introduced 
together with the contact interaction and require the introduction of another isoscalar term. The 
BCS-BEC crossover of neutrons pairs in symmetric and asymmetric nuclear matter is studied by 
using these contact interactions. It is shown that the bare and screened pairing interactions lead to 
different features of the BCS-BEC crossover in symmetric nuclear matter. For the screened pairing 
interaction, a two-neutron BEC state is formed in symmetric matter at kpn ~ 0.2 fm~^ (neutron 
density p„/po ^ 10~^). Contrary the bare interaction does not form the BEC state at any neutron 
density. 
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I. INTRODUCTION 



Low density nuclear matter tends to be much more 
interesting than a simple zero density limit of the bulk 
physics of nuclei. Indeed, new phenomena show up and 
are mainly associated with the formation of bound states 
or with the emergence of strong correlations [l|, [3, 0, 

Q. For instance, in symmetric nuclear matter, neu- 
trons and protons become strongly correlated while the 
density decreases and a deuteron BEC state appears at 
very low density Q ■ The deuteron- type correlations give 
extra binding to the nuclear equation of state and induce 
new features at low density [8|. This transition belongs 
to the BCS-BEC crossover phenomena which have been 
extensively studied in several domains of physics and has 
recently been experimentally accessible in cold atomic 
gas (see Ref. @ and references therein). 

In nuclear matter, neutron pairs are also strongly cor- 
related. Theoretical predictions suggest that, at density 
around po/10 where po=0-16 fm~'^, the pairing gap 
may take a considerabl y la rger value than that at nor- 
mal nuclear density po The density dependence of 
the pairing gap at low density is unfortunately not yet 
completely clarified and still awaits a satisfactory solu- 
tion [ll|, Therefore, it could be interesting to 
explore pairing interactions based on Brueckner theory 
and its consequences to the BCS-BEC crossover. Indeed, 
pairing at low density is relevant for different purposes: 
for the understanding of neutron-rich exotic nuclei near 
the drip line [11 [S [lE [H [S [11 , where the long tails 
of density profiles give rise to "halo" or "skin" behavior, 
or for the expanding nuclear matter in heavy ion colli- 
sions or even for the physics of neutron stars, where 
several physical phenomena, such as cooling and glitches, 



are thought to depend very sensitively on the size and 
the density dependence of the pairing gap [2l|, [13, IH • In 
^^Li, the wave function of the two neutrons participating 
to the halo nucleus has been analyzed with respect to the 
BCS-BEC crossover [2J|. It has been shown that as the 



distance between the center of mass of the two neutrons 
and the core increases, the wave function changes from 
the weak coupling BCS regime to the strongly correlated 
BEC regime. This is due to the fact that the pairing 
correlations are strongly density-dependent [l^ and the 
distance between the two neutrons and the core provides 
a measure of the pairing strength. 

It should be emphasized that the bare nuclear interac- 
tion in the particle-particle channel should be corrected 
by the medium polarization effects [13, [ll| (usually re- 
ferred to as the screening effects). These effects have been 
neglected for a long time since the nuclear interaction is 
already attractive in the ^Sq channel without the medium 
polarization effects, contrary to the Coulomb interaction 
for which the medium polarization effects are absolutely 
necessary to get an attractive interaction between elec- 
trons. However, several many-body methods have been 
developed recently to include the medium polarization 
effects in the calculation of the pairing gap such as a 
group renormalization method [20|, Monte-Carlo meth- 
ods calculation [27,[2l, [11 and extensions of the Brueck- 
ner theory (Tol. [ll | . These calculations, except the one 
presented in Ref. 27[ , predict a reduction of the pairing 
gap in neutron matter. 

Note that, based on the nuclear field model, it has also 
been suggested that the medium polarization effects con- 
tribute to the pairing interaction in finite nuclei and in- 
crease the pairing gap [30, 31, 32]. To understand this ap- 
parent contradiction between neutron matter and finite 
nuclei, a Brueckner calculation including the medium po- 
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larization effects in both symmetric and neutron matter 
has been performed in Ref. It has been shown that 
the medium polarization effects are different in neutron 
matter and in symmetric matter. The medium polariza- 
tion effects do not reduce the pairing gap in symmetric 
matter, contrary to that in neutron matter. Instead, in 
symmetric matter, the neutron pairing gap is much en- 
larged at low density compared to that of the bare calcu- 
lation. This enhancement takes place especially for neu- 
tron Fermi momenta kp^ < 0.7 fm~^. This could explain 
why the medium polarization effects increase largely the 
pairing correlations in finite nuclei but decrease it in neu- 
tron matter. 

In this paper, we propose an effective density- 
dependent pairing interaction which reproduces both the 
neutron-neutron (nn) scattering length at zero density 
and the neutron pairing gap in uniform matter obtained 
by a microscopic treatment based on the nucleon-nucleon 



interaction pj|. The proposed interaction has isoscalar 
and isovector terms which could simultaneously describe 
the density dependence of the neutron pairing gap for 
both symmetric and neutron matter. Furthermore, we 
invent different density-dependent interactions to de- 
scribe the "bare" and "screened" pairing gaps, together 
with the asymmetry of uniform matter, given in Ref. [ll| . 
Then, we explore the BCS-BEC crossover phenomena in 
symmetric and asymmetric nuclear matter. 

This paper is organized as follows. In Sec. [TTl we dis- 
cuss how to determine the isoscalar and isovector density- 
dependent contact interactions. Applications of those in- 
teractions to the BCS-BEC crossover are presented in 
Sec. mil We give our conclusions in Sec. IIVI 



II. DENSITY-DEPENDENT PAIRING 
INTERACTION 

Recently, spatial structure of neutron Cooper pair 
in low density nuclear matter has been studied using 
both finite range interactions like Gogny or GSRS and 
density-dependent contact interactions properly adjusted 
to mimic the pairing gap obtained with the former inter- 
actions [2^. It was found that the contact interactions 
provide almost equivalent results compared to the finite 
range ones for many properties of the Cooper pair wave 
functions. It is thus reasonable to investigate the evolu- 
tion of the Cooper pair wave function with respect to the 
density and the isospin asymmetry using contact inter- 
actions adjusted to realistic interactions. In this paper, 
we take a contact interaction w„„ acting on the singlet 
^S channel, 



\k') = 



1 



-VQg[Pn,Pp]d{k,k') 



(1) 



where the cutoff function 9{k, k') is introduced to remove 
the ultraviolet divergence in the particle-particle channel. 
A simple regularization could be done by introducing a 
cutoff momentum kc- That is, 0{k,k') = 1 if fc, fc' < kc 



and otherwise. In finite systems, a cutoff energy Cc is 
usually introduced instead of a cutoff momentum kc- The 
relation between the cutoff energy and the cutoff momen- 
tum may depend on the physical problem, and it is known 
that the pairing strength vq depends strongly on the cut- 
off. A detailed discussion on the different prescriptions 
used in the literature are then presented in Appendix [Al 
In this paper, we choose the prescription 3 in the Ap- 
pendixl^so that the adjusted interaction can be directly 
applied to Hartree-Fock-Bogoliubov calculations. 

In Eq. Ill, the interaction strength vq is determined 
by the low energy scattering phase-shift, that fixes the 
relation between vq and the cutoff energy Cc, while the 
density-dependent term g[pn, Pp] is deduced from predic- 
tions of the pairing gaps in symmetric and neutron mat- 
ter based on the nucleon-nucleon interaction [11.] . The 
density-dependent term accounts for the medium effects 
and satisfies the boundary condition g 1 for p —f 0. 
The volume type and surface type pairing interactions 
have g — 1 and g — at p — pQ, respectively. In this 
paper, we introduce more general types of pairing inter- 
actions and the novelty is a dependence on the ratio of 
neutron to proton composition of the considered system. 
We thus define a function 



gi[Pn,Pp] = 1 - fs{I)Vs 



PO 



- fn{I)nn 



,(2) 



where / is the asymmetry parameter, defined as / = 
{N ~ Z)/{N + Z), and po=0.16 frn'^ is the saturation 
density of symmetric nuclear matter. We insert the func- 
tion gi into Eq. ([T]) as 17 = (?i . The goal of the functional 
form in Eq. ^ is to reproduce the theoretical calcula- 
tion of the pairing gap in both symmetric and neutron 
matter and also to be used for prediction of the pairing 
gap in asymmetric matter. It could also be applied to 
describe pairing correlations in finite nuclei by acquiring 
an explicit dependence on the coordinate r in the den- 
sity p{r) and the asymmetry parameter I{r). In Eq. ([2]), 
the interpolation functions and fnil) are not ex- 

plicitly known but should satisfy the following condition 
/,(0) = /„(1) = 1 and fs{l) = /„(0) = 0. The density- 
dependent function gi is flexible enough and we can ob- 
tain an effective pairing interaction which reproduces the 
density dependence of the pairing gap in uniform mat- 
ter. It should however be noticed that the interpolation 
functions and fnil) cannot be deduced from the 

adjustment of the pairing gap in symmetric and neutron 
matter. For that, theoretical calculations in asymmetric 
nuclear matter or application to exotic nuclei are neces- 
sary. In this paper, we choose fs{I) = 1 — fn{I) and 
/„(/) = /. Many different interpolation functions could 
be explored but we think that it has little consequences 
to the BCS-BEC crossover. 

We have also explored other density-dependent func- 
tionals, introducing an explicit dependence on the isovec- 
tor density, I — r/s (^-^^ — Vi '^"p^^'' ) ' introduc- 



ing the neutron and proton densities, 1 — ?]„ f ^ 
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FIG. 1: (Color online) Phase shifts for s-wave nucleon- 
nucleon scattering as a function of the center of mass energy. 
In the left panel are shown nn phase shifts obtained from 
Argonne i;i8 potential [3l (stars) and the result of the best 
adjustment obtained with a contact interaction for a set of 
cutoff energies (from 10 to 120 MeV). In the right panel are 
shown the s-wave phase shifts for various channels: nn, np 
and pp. The np and pp phase shifts have been provided by 
the Nijmegen group (http: / / nn-online.orgp . 



T7p • The isospin dependence of those functionals 

is fixed but these functional forms are not flexible enough 
so that the adjustment becomes sometimes impossible. 
For instance, the pairing gap in symmetric and neutron 
matter including the medium polarization effects and cal- 
culated in Rcf. [ill ] cannot be reproduced with such func- 
tionals. In the following, we have therefore considered 
only the density-dependent interaction given in Eq. 

Recently, another attempt has been done to introduce 
an isospin dependent volume type pairing interaction 
in a very different approach [34| . This interaction has 
been adjusted to reproduce empirical mass formula over 
a thousands of nuclei, but the pairing gap calculated with 
this pairing interaction compares very badly with realis- 
tic calculations in uniform matter. 



A. The free interaction 

In Ref. 14], it has been proposed to deduce the free 
interaction parameter vq from the low energy phase shift 
of nucleon-nucleon scattering. The nn, np and pp phase 
shifts versus the center of mass energy are shown on the 
right panel of Fig.[TJ It is clear from this figure that each 
of these three channels are different and the interaction 
parameter vq should depend on the channel of interest. 
In this paper, we are interested only in the nn channel. 
We then express the phase shift as a function of the cutoff 
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TABLE I: For a given cutoff energy ec, the parameters r„„, 
a and vo are determined by the scattering length arm which 
reproduces the phase shift in the low energy region ec.m. < 
2 MeV. The interaction strengths Vq and Vq" are obtained 
from the empirical value of a„„=-18.5 fm and from the unitary 
limit, defined as dnn ^ OO. 



momentum kc and the scattering length o„„ p^ . [s^ ] : 



fccot(5 



QfTT 

1 



1 -I- akc 



ak , kr ~ k 
— In-^ 

2 kc + k 



In — r 



(3) 
(4) 



where a = 2o„„/(7r — 2fccO„„). In this way, for a given 
cutoff momentum kc, the phase shift can be adjusted us- 
ing the scattering length a„„ as a variable. The results 
are shown in Fig. [T] and Table U for a set of cutoff en- 
ergies Cc — h'^kl/m (note that we use the reduced mass 
m/2). We have found that, for cutoff energies larger than 
20 MeV, the optimal parameters cannot reproduce the nn 
phase shift in the low energy region [ec.m. <2 MeV) and 
in the higher energy region simultaneously. We mainly 
focus on the adjustment of the nn phase shift in the 
low energy region. At higher energies, or equivalently 
at higher densities in nuclear matter, the medium effects 
modify the interaction anyhow and generate a density 
dependent term in the interaction. 

Fixing Cc and a„n, one can deduce the effective range 
fnn — 4/7rfcc, the parameter a and the interaction 
strength vq = ^li^af? jm. The values of those param- 
eters are given in Table HI The value of the free inter- 
action parameter deduced from the empirical value 
of the scattering length a„„ = —18.5 fm is also indi- 
cated. One can see that the difference between vq and 
fg is small, as much as 3%. Indeed, as we are in a 
regime of large scattering length, one can deduce the in- 
teraction strength approximately from the relation wg ~ 
v"^ (1 H- -K/2kcann + • ■ • ) where — ~2n'^tf/mkc is the 



interacting strength in the unitary limit {kca„ 



oo). 



B. The density-dependent function g[pn,pp] 

The density-dependent function g is adjusted to repro- 
duce the pairing gaps in symmetric and neutron matter 
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obtained from Ref . flT'l . Pairing in uniform nuclear mat- 
ter is evaluated with the BCS ansatz: 



\BCS) 



(5) 



fe>0 



where Uk and Vk represent the BCS variational pa- 
rameters and a|,| are creation operators of a particle 
with momentum k and spin "f on top of the vacuum 
I-) [35, 36, 37]. The BCS equations are deduced from 
the minimization of the energy with respect to the vari- 
ational parameters Uk and Vk- For a contact interaction, 
the equation for the pairing gap A„ takes the following 
simple form at zero temperature, 



VQ9[Pn,Pp] 

2(27r)3 



En{k) 



e{k,k) , 



(6) 



where 9{k, k) is the cutoff function associated to the con- 
tact interaction £'„(fc) = ^ (e„(fc) — Vn)'^ + is the 
neutron quasi-particle energy, e„(fc) = ft^fc^/2m* is the 
neutron single particle kinetic energy with the effective 
mass m*. We define the effective neutron chemical po- 
tential — iin~Un, where the neutron mean field poten- 
tial Un is subtracted from the neutron chemical potential 
Hn- The effective neutron chemical potential Vn gives the 
neutron density, 



2 
V 



(7) 



where V is the volume and nn{k) is the occupation prob- 
ability defined as 



1 - 



En{k) 



(8) 



Finally, the neutron Fermi momentum kpn is defined as 
p„ = 4„/37r2. 

We have chosen to adjust our interaction to the re- 
sults of nuclear matter pairing gaps in Ref. since 
it is the only calculations performed for both symmet- 
ric and neutron matter. We adjust the contact pair- 
ing interaction so that it reproduces the position and 
the absolute values of the maxima of the pairing gaps 
in symmetric and neutron matter. For the bare pairing 
gap, the maximum is located at kpn — 0.87 fm"-'^ with 
A„=3.1 MeV for both symmetric and neutron matter, 
while for the screened pairing gap, the maximum is at 
kpn ~ 0.60 fm~^ with A„=2.70 MeV for symmetric mat- 
ter and kpn — 0.83 fm"-'^ and A„=1.76 MeV for neutron 
matter. The value of the parameters t]s and 7y„ are freely 
explored in the real axis while the parameters ols and 
an are imposed to be positive to avoid singularities. The 
neutron effective mass to* is obtained from SLy4 Skyrme 
interaction since it is widely used in nuclear mean-field 
calculations. The results of the fits are given in Table HIl 
and the pairing gaps are shown in Fig. [51 

One should note that for the bare interaction, even 
if the pairing gap is identical in symmetric and neutron 
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bare 


60 MeV 


0.598 
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0.525 




20 MeV 
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10 MeV 
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60 MeV 
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40 MeV 
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10 MeV 
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6U MeV 
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n oo 
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5 = 5i + 52 


40 MeV 


1.80 


0.27 


1.61 


0.122 


7?2 = 0.8 


20 MeV 


2.06 


0.31 


1.70 


0.122 




10 MeV 


2.44 


0.37 


1.66 


0.0939 



TABLE II: Parameters of the function g defined in Eq. 
These parameters are obtained from the fit to the pairing 
gaps in symmetric and neutron matter. These are the param- 
eters obtained from the adjustment of the bare gap and the 
screened gap with g = gi, and the screened gap including the 
additional function g2. The efi'ective mass is obtained from 
SLy4 Skyrme interaction. Note that Ec is the cutoff for the 
quasi-particle gap equation ® while ec is that for the two- 
body scattering so that Ec — ec/2. See the text for details. 



matter, the adjusted contact interaction is not necessar- 
ily isoscalar. Indeed, the transformation from the Fermi 
momentum, the x-axis of Fig. [51 to the density is dif- 
ferent in symmetric nuclear matter, p/po = ikpn/kpo)^ 
(where po = 2/(37r^)fc|,Q=0.16 fm~'^), and in neutron 
matter, p/po = 0.5{kFn/kFo)'^ ■ Therefore, an interac- 
tion which depends only on the ratio p/ po gives different 
results if it is plotted as a function of kpn in symmet- 
ric and neutron matter. As the pairing gap calculated 
with the bare interaction pd| is quasi-identical in sym- 
metric and neutron matter when it is plotted versus k pn , 
one can then deduce the following relations between the 
parameters of the density-dependent term gi (neglecting 
the isospin dependence of the effective mass): ag — an 
and ris = 77„/2"". 

For the bare pairing gap and for a given cutoff energy 
Ec, the position and the maximum value of the gap are re- 
produced well by the contact interaction in Fig. [51 How- 
ever, in the high Fermi momentum region kpn > 1 fm~^, 
we can see appreciable difference between the microscopic 
predictions and the pairing gap obtained from the con- 
tact interactions. The best agreement is obtained for a 
cutoff energy Ec = 40 MeV. In the screened case, the de- 
pendence of the pairing gap on kpn is badly reproduced, 
especially for symmetric nuclear matter. This is because 
the maximum position of the pairing gap is shifted to- 
wards a lower neutron Fermi momentum (one third in 
density from that for the bare gap). Consequently, the 
density dependence of the function gi becomes stiffer in 
the "screened" case than in the bare case, and the gap 
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FIG. 2; (Color online) Pairing gap in symmetric, asymmetric 
{xp — pp/p) and neutron matter adjusted to the "bare gap" 
(upper panel) or to the "screened gap" (lower panel) for var- 
ious cutoflt energies Ec. The pairing gap calculated from the 
microscopic treatment presented in Ref. is also shown as 
star symbols. 



drops faster after the maximum, as it is shown in Fig. [2l 
This may indicate that the screened interaction has a 
different density dependence and cannot be cast into a 
simple power law of the density as in Eq. Indeed, 
in Ref. jTlj, the medium polarization effects have been 
analyzed at the level of the interacting potential, and it 
was shown that the medium polarization effects emerge 
at very low density and remain relatively constant. To 
simulate such effects, it seems necessary to introduce a 
new term, 52 in Eq. ([2|). We propose for 32 a simple 
isoscalar constant which switches on at a very low value 
of the density {kp ~ 0.15 fm^-'^) and switch off around 
the saturation density. The following form satisfies this 
condition, 



92 = m 



( 1 + e 05 1 



- 1 



(9) 



The new pairing interaction with g — gi + 92, hereafter 
named screened-II, has then only one new adjustable pa- 
rameter 772. As the medium polarization effects could 
also change the density-dependent term gi, all the 5 pa- 
rameters have to be re-adjusted. In Tables |TT] and IIIIl 
we give the new parameters t^s, c^s, Vn and «„, obtained 
for several values of 772. The cutoff energy is fixed to 
be i?c =40 MeV. The corresponding pairing gap is repre- 
sented in Fig. [21 The best fit is obtained for the value 
7^2=0.8. Eq. ^ may not be a unique way to take into 
account the medium polarization effects. Nevertheless it 
is simple enough to apply to the BCS-BEC crossover, so 
this is why we adopt this functional form. 

Solving the gap Eq. the neutron effective chemi- 
cal potential Vn is determined for a given interaction at 




FIG. 3: (Color online) Pairing gap calculated in symmetric, 
asymmetric and neutron matter with the screened-II interac- 
tion {g — <7i + 32) and for several values of 772 as indicated 
in the legend. The corresponding parameters are given in 
Table Hill See the text for details. 



a given neutron Fermi momentum kpn- The neutron ef- 
fective mass m* , the effective neutron chemical potential 
i^n — ^J'n — Un and the difference Vn — ^Fn are represented 
in Fig. m as a function of the neutron Fermi momentum 
in symmetric and neutron matter. The neutron effective 
mass and the neutron potential Un are deduced from the 
SLy4 Skyrme interaction. Note that the neutron density 
Pn is changed into the neutron Fermi momentum, kpn- 
We have selected from Table ITTl the bare and the screened- 
II pairing interactions for a cutoff energy Ec = 40 MeV. 
In the absence of pairing, the effective neutron chemi- 
cal potential Vn is identical to the neutron Fermi kinetic 
energy, Vn = ipn where epn = in{kpn)- The difference 
Vn—^Fn is then null in the absence of pairing correlations, 
otherwise it is negative as shown in Fig.|4l From this dif- 
ference, one can estimate the relative importance of the 
pairing correlations: in neutron matter the screened-II 
interaction leads to weaker pairing correlations compared 
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TABLE III: Parameters of the the screened-II interaction, 
where the density-dependent function of the pairing interac- 
tion is taken tohe g = g\+g2. The functional forms g\ and g2 
are obtained by fitting the screened pairing gap for several val- 
ues of 772. The energy cutoff is taken to be Ec = 40 MeV and 
the neutron effective mass is deduced from the SLy4 Skyrme 
interaction. 
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FIG. 4: (Color online) Comparison of the neutron effective 
mass m'^, the effective chemical potential iy„ = fi„ — Un, 
and the difference iy„ — epn calculated with the bare and the 
screened-II contact interactions. The parameters of the pair- 
ing interactions are taken from Table lUwith the cutoff energy 
Ec=40 MeV. The effective mass mj^ and the neutron potential 
U„ are taken from SLy4 Skyrme interaction. The two arrows 
indicate the lower and upper limits for the condition v„ < 
with the screened-II interaction in symmetric nuclear matter. 



to the bare one, while in symmetric matter, the screened- 
II interaction give much stronger pairing correlations for 
kpn < 0.7 fm~-^ and less for kpn > 0.7 fm^^. 

It is easy to show that the gap © and the occupation 
probability ^ go over into the Schrodinger-like equation 
for the neutron pair wave function "if pair [38[ . 

p2 1 

~'^pair + [1 - 2n„(fc)]— Tr Vnn'^pair = 2Vn'^pair- (10) 

See Eq.( [T3|) in Sec. IIIII for a proper definition of the 
neutron pair wave function ^pair- Notice that, at zero 
density where nn{k) — 0, Eq. pO|) is nothing but 
the Schrodinger equation for the neutron pair. From 
Eq. pU|) . one usually relates the effective neutron chemi- 
cal potential zy„ to be a half of the "binding energy" of a 
Cooper pair. The Cooper pairs may then be considered 
to be strongly correlated if i^n is negative. Notice from 
Fig. m that the effective neutron chemical potential i/n 
becomes negative with the screened-II interaction in sym- 
metric nuclear matter for fci?„=0.05-0.35 fm^^, but not 
at all in neutron matter. One could then expect in this 
case that the Cooper pair may resemble a closely bound 
system (BEC) in symmetric nuclear matter at very low 
density while in neutron matter it should behave like that 
of the weak coupling BCS regime. However, to go beyond 
this rough interpretation, we need to study more accu- 
rately the BCS-BEC crossover in asymmetric matter. 
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-0.77 



BCS boundary 
unitarity limit 
BEC boundary 



TABLE IV: Reference values of {k 

An/epn and Vn/eFn characterizing the BCS-BEC crossover in 
the regularized model for the contact interaction. The values 
dn, P{d„), and ^mis are the average distance between neu- 
trons dn — pn^^'' , the probability for the partner neutrons 
correlated within the relative distance d„, and the rms radius 
of Cooper pair , respectively. The numbers have been taken 
from Refs. [251. l4g|. See the text for details. 



III. APPLICATION TO THE BCS-BEC 
CROSSOVER 

Going from the weak coupling BCS regime, around the 
saturation density po, down to the BCS-BEC crossover, 
for densities < po/lQ, it has been shown that the spa - 
tial structure of the neutron Cooper pair changes [25j . 
It is indeed expected that the correlations between two 
neutrons get large as the density decreases and as a con- 
sequence, the BCS-BEC crossover occurs in the uniform 
matter at low density. However, being of the second or- 
der, this transition is smooth. Hereafter, we clarify the 
boundaries of the BCS-BEC phase transition by using a 
regularized gap equation. 

Although the BCS ansatz ^ has been developed to 
describe the Cooper pair formation in the weak BCS 
regime [35*1, it has been shown that the BCS equations are 
also valid in the strong BEC condensation regime [38:,, 3^. 
The BCS equations are thus adopted as a useful frame- 
work to describe the intermediate BCS-BEC crossover 
regime at zero temperature [i^l • It has been proposed to 
define the limit of the BCS-BEC phase transition using a 
regularized model for the pairing gap [1^, HO, Elj . In this 
model, the BCS gap (O is combined with the relation be- 
tween the interaction strength and the scattering length 
which has a similar divergent behavior. The difference 
between those two divergent integrals gives a regularized 
equation. 



47ra„„ 2V \En{k) en{k) 




which has no divergence. The gap equation ([6]) can be 
solved analytically for the contact interaction with a con- 
straint of the particle number conservation ([7]) . The solu- 
tion of this regularized gap equation is independent of the 
strength of the interaction, and the gap is uniquely de- 
termined by the value of the scattering length a„„. From 
Eq. pT|) . one can study the boundaries of the BCS-BEC 
phase transition with respect to the dimensionless order 
parameter kpnCinn- We give in Table HVl the values of 
several quantities which specify the phase transition: the 
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probability P{dn) for the partner neutrons to be corre- 
lated within the relative distance (i„ (d„ is the average 

— 1/3 

distance between neutrons dn = Pn ) , the ratio of the 
rms radius to the mean neutron distance £,rms/dn, the ra- 
tio of the pairing gap to the single particle kinetic energy 
^n/f-Fn and also the ratio of the effective neutron chemi- 
cal potential to the single particle kinetic energy Vn/^Fn- 
As we already mentioned, these boundaries are indicative 
because the phase transition is smooth at the boundaries 
being of the second order. For instance, even if the nu- 
clear matter does not enter into the BEC regime, we will 
show that the Cooper pair wave function is already very 
similar to the BEC one when it is close. 

A drawback of this regularized model is that the rela- 
tion between the dimensionless order parameter kFnO-nn 
and the density of the medium is unknown. To relate the 
order parameter to the density, one has to re-introduce 
the pairing strength in the gap equation ©. We could 
consider for instance a contact interaction with a cutoff 
regularization. The density will then trigger the phase 
transition for a given pairing interaction strength. In 
the following, we study the BCS-BEC phase diagram in 
asymmetric nuclear matter for the two pairing interac- 
tions discussed in Sec. [TTl Namely we explore the prop- 
erties of the Cooper pair wave function obtained by the 
bare and the screened-II interactions presented in Ta- 
ble |TT] for a fixed cutoff energy i?c=40 MeV. 

The BCS approximation provides the Cooper pair 
wave function ^pair{k) [ssl. [Sa. [37j 

■^pa^r{k) = C {BC S\a) {k X)d\-k [)\BC S) (12) 
= CukVk ■ (13) 

The radial shape of the Cooper pair wave function is de- 
duced from the Fourier transform of u^Vk = A„/2£'„(A:) 
in Eq. (fT^. The rms radius of Cooper pairs is then 

given by ^„ns = /PI = J drr'^\-i'pa,r{rW ■ The 
rms radius ^rms and the Pippard's coherence length, 
— h'^kpn/m'^TTAn, give similar size of the Cooper pair 
in the weak coupling regime. The rms radius ^rms is nev- 
ertheless a more appropriate quantity in the BCS-BEC 
crossover region as well as in the strong BEC coupling 
region. In order to estimate the size of Cooper pairs a 
reference scale is given by the average distance between 
neutrons c?„. If the rms radius of Cooper pairs is larger 
than dn, the pair is interpreted as an extended BCS pair 
while the Cooper pair will be considered as a compact 
BEC pair if the rms radius is smaller than the average 
distance. Let us introduce another important quantity 
which also gives a measure of the spatial correlations: the 
probability P(r) for the partners of the neutron Cooper 
pair to come close to each other within the relative dis- 
tance r. 



P(r) = / dr'r'■'\^^Jpa^r{r')\' 



(14) 



The order parameters listed in Table IIVI are closely re- 
lated. For instance, approximating ^rms by ^p, it could 



easily be shown that the ratio ^rms/dn is proportional to 
epn/ An- Then, the strong coupling regime is reached if 
the ratio An/epn is large. The order parameter t'„, the 
effective neutron chemical potential, could be interpreted 
as a half of the binding energy of Cooper pairs at finite 
density according to the Schrodinger-like Eq. (fTUl) . 

As shown in the Appendix[B]it is convenient to decom- 
pose the Cooper pair wave function into 



'i'pairir) = «'i(r) + ^'2(r) 



where 



*2(r) = C 



„ fc2 sin(fcr) 
dk — -, 

En{k) kr 

' fc^ sin(fcr) 
En{k) kr 



(15) 

(16) 
(17) 



Choosing k^o/k^ 3> 1, with fco — -\/2to* j/n/^i, it is possi- 
ble to find an analytic expression for '52: 



*2(r) 



C 2m* 
r h? 



si(fcoor) 



(18) 



where si(u) is the sinus integral defined as si(u) = 
J°° dz [sin(z)/2]. It is clear from Eq. that the term 
^2 has a 1/r-type singularity. This singularity is due 
to the nature of the contact interaction which does not 
contain a hard core repulsion. With the hard core repul- 
sion, the wave function goes to zero as r — + [2^. In 
the outer region (r > 3 fm), the wave function behaves 
in the same way if the contact interaction is deduced 
properly from the microscopic calculations. We checked 
the convergence of the wave function with respect 
to the parameter k^o- We found that the convergence is 
reached with fcoo ~ 2fcn as is is shown in Fig. [TT] in the 
Appendix [B] In Ref. [23, Matsuo introduced the cutoff 
momentum kc to calculate the pair wave function (jlSp . 
We have compared the pair wave function "^pairif) with 
the one obtained by Matsuo. The two wave functions 
give essentially the same results, except for the low den- 
sity region. In the worst case, the wave function of Mat- 
suo's treatment increases the rms radius by about 10% 
compared to the one obtained by the wave function (fTSjl . 

The neutron Cooper pair wave function r'^\'^pair{f)\'^ 
is shown in Fig. [5] as a function of the relative distance 
r between the pair partners taking different Fermi mo- 
menta fcFn=l-l, 0.8, 0.5 and 0.2 fm~^, which correspond 
respectively to the densities: pn/po—0.3, 0.1, 0.03 and 
0.002. Calculations in symmetric, asymmetric and neu- 
tron matter are shown in the left, middle and right pan- 
els, respectively. In Fig.[5l we observe that the spatial ex- 
tension and the profile of the Cooper pair varies strongly 
with the density. A large extension is found close to the 
saturation density at fci?„ = l.l fm^^. The profile of the 
wave function behaves as an oscillation convoluted by a 
decreasing exponent and casts into the well known limit 
^ Ko{r/TT£^p)sm{kFr)/kFr 35|. This indicates that the 
Cooper pair is in the weak coupling BCS regime. At lower 
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symm. matt. asymm. matt. (x=0.3) 



neut. matt. 




5 10 15 20 5 10 15 20 5 10 15 20 

r [fm] r [fm] r [fm] 



FIG. 5: (Color online) Neutron Cooper pair wave function r^\'i/pair{r)f as a function of the relative distance r between the 
pair partner at the Fermi momenta fcFn = l-l, 0-8, 0.5 and 0.2 fm~^. 



densities, the Cooper pair shrinks and the osciUation dis- 
appears. The wave function resembles now the strong 
couphnghmit (BEC) - exp{- y/ im/ \fi\r)/r This 
is an indication that a possible BCS-BEC crossover may 
occur in uniform matter. 

It should be remarked that the latter limit seems well 
pronounced in symmetric matter with the screened gap 
(see the panel at the left bottom corner of Fig. [5]). We 
show in Fig. [6] the evolution of the occupation proba- 
bility ([8]) in symmetric matter for the two pairing in- 
teractions. For the screened-II interaction, the pairing 
correlations becomes strong at low densities as the occu- 
pation probability is considerably different from the step 
function. In the case of the bare interaction, the correla- 
tions are not so strong to change n„(fc) drastically even 
at low densities. It should be noticed that this analysis is 
independent of the detailed structure of the Cooper pair 
wave function. This change of the occupation probability 
proves that the behavior of the Cooper pair wave function 
is not an artifact induced by the zero range behavior of 
the contact interaction but indeed is physical. It is clear 
that low density symmetric nuclear matter is much more 
correlated with the screened-II interaction than with the 
bare one. This is also the case for the BCS-BEC crossover 
as will be discussed below. 

Let us now discuss the BCS-BEC crossover which may 
depend on the pairing interactions and also on the asym- 



metry of the nuclear medium. In the following, we study 
the different order parameters in Table HVl for the bound- 
aries of the BCS-BEC phase transition. Fig. [7] shows the 



bare interaction 



screened-II interaction 




■■'■)•'..' ' .1 

k =0.8-fm 

Fn 


k„ =1.1 fm 

Fn 


- k„ =0.5 fm"' ■■ 




\ 
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I 


=0.2 fm'' \ 

Fn 

■ — 




. \ - 
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, •,■ . 



0.5 1 

e„(k)/e„„ 



0.5 1 



FIG. 6: (Color online) Occupation probability n„(fc) for sym- 
metric matter defined by Eq. ([8} as a function of the ra- 
tio of the single particle kinetic energy to the Fermi energy 
en(fc)/eFn for a set of Fermi momenta kFn= 1-1, 0.8, 0.5 and 
0.2 fm~^. We compare the results of the bare interaction (left 
panel) with the screened-II interaction (right panel). The 
values of the effective neutron chemical potential Vn/epn are 
indicated by the filled circles. 
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FIG. 7: (Color online) Probability P for the partner neutrons within two typical lengths, 3 fm and ci„, as a function of the 
neutron Fermi momentum kpn in symmetric (left panel), asymmetric (central panel) and neutron matter (right panel). The 
boundary of the BCS-BEC crossover is denoted by the dashed line. 



probabilities P{r) for the partner neutrons to be corre- 
lated within the typical scales, r=3 fm and r—dn- The 
former scale is the typical range of the nucleon-nucleon 
force. For the bare interaction, the probability P(3 fm) 
has similar behavior in symmetric and asymmetric mat- 
ter as a function of kpn. For the screened-II interaction, 
there is a noticeable isospin dependence. A low den- 
sity shoulder appears in symmetric matter, at around 
kpn ~ 0.25 fm~^ [Pn/ Pq ^ 0.003). Then, it becomes 
smaller as the asymmetry increases and eventually dis- 
appears in neutron matter. In neutron matter, the strong 
concentration of the pair wave function within the inter- 
action range 3 fm, P(3 fm) > 0.5, is realized in the den- 
sity region kpn ^ 0.3—1.1 fm~^ (or Pn/po ~ 0.007 — 0.3) 
for both pairing interactions. For symmetric matter, on 
the other hand, this region is different for the two pairing 
interactions: the strong correlation occurs at much lower 
density region for the screened-II interaction than for the 
bare one. This property can also be confirmed by the 
probability P{dn). For the two pairing interactions, the 
Cooper pairs in symmetric and asymmetric matter enter 
into the crossover regime at almost the same density. The 
crossover in neutron matter occurs somewhat at lower 
density for the screened-II interaction. As the density 
decreases, a different behavior is observed between the 
two pairing interactions for symmetric matter. While the 
probability P{dn) decreases and goes back to the weak 
BCS regime for the bare interaction at very small density 



below kpn ~ 0.1 fm""'^, the probability P{dn) continue to 
increase up to 1 for the screened-II interaction at very low 
densities, kpn < 0.7 fm~^ {p — n/p = < 0.07 fm~^). 

We study further the BCS-BEC crossover by looking 
at the rms radius ^rms and the neutron pairing gap A„. 
In Fig. [51 we show the rms radius ^rms as a function of 
the neutron Fermi momentum kpn as well as the order 
parameter ^rms/dn- The rms radius of the Cooper pair 
is less than 5 fm in the region kpn ~ (0.4 — 0.9) fm~^ 
{Pn/po ^ 0.01 — 0.15) in the three panels for the bare 
interaction. The screened-II interaction gives different 
effects in symmetric and asymmetric matter: it increases 
the rms radius for the neutron matter, while the rms ra- 
dius stays small around 4 fm even at very low density at 
kpn ~ 0.15 fm^^ (Pn/Po ^ 0.0007) in symmetric mat- 
ter. In the lower panels is shown the ratio of the rms 
radius to the average distance between neutrons d„. For 
the bare interaction, the size of the Cooper pair becomes 
smaller than d„ for the Fermi momentum kpn < 0.8 fm"-'^ 
(pn/po ~ 0.1) in general. There are substantial differ- 
ences for symmetric and asymmetric matter in the case 
of the screened-II interaction. The crossover region be- 
comes smaller for the neutron matter, while the crossover 
region increases in the cases of asymmetric (xp — 0.3) 
and symmetric matter. Especially, the correlations be- 
comes strong in symmetric matter and the Cooper pair 
reaches almost the BEC boundary at kpn ^ 0.2 fm~^ 
{Pn/po ^ 0.002). Notice that the two neutrons system is 
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FIG. 8: (Color online) Top panels: Comparison between the rms radius ^rms of the neutron pair and the average inter-neutron 
distance dn — p^^^'^ (thin line) as a function of the neutron Fermi momentum kpn in symmetric (left panel), asymmetric 
(central panel) and neutron matter (right panel). Bottom panels: The order parameter ^rms/dn as a function of kpn- The 
boundaries of the BCS-BEC crossover are represented by the two dashed lines, while the unitary limit is shown by the dotted 
line. The two pairing interactions are used for the calculations. 




FIG. 9: (Color online) Ratios A„/eFn and v„/^Fn plotted as a function of the neutron Fermi momentum kpn in symmetric (left 
panel), asymmetric (central panel) and neutron matter (right panel). The boundaries of the BCS-BEC crossover are shown by 
the two dashed lines, while the unitary limit is given by the dotted line. See the text for details. 
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known experimentally to have a virtual state in the zero 
density limit. We have shown that this virtual state could 
lead to a strongly correlated BEC state at low density in 
symmetric nuclear matter according to the screened-II 
interaction. 

The two other order parameters l^n/^Fn and VnjtFn 
are shown in Fig. [9l These results confirm the BCS-BEC 
crossover behavior which is found in Fig. [51 Namely, 
in symmetric matter, the gap A„ is much enhanced by 
the screened interaction in the low density region, while 
no enhancement can be seen in neutron matter. As ex- 
pected, the effective chemical potential Vn induced by 
the screened-II interaction becomes negative for kpn ^ 
0.05-0.3 fm-i [pn/po 0.00002-0.01) in symmetric nu- 
clear matter. This strong correlation is reduced in asym- 
metric matter and is almost absent in the neutron matter 
as can be seen in Fig. [HI It should also be remarked that 
the order parameter VnjtFm due to the Schrodinger-like 
Eq. pO)) . gives the same BCS-BEC crossover behavior as 
that of the order parameter ^rms/dn- The neutron effec- 
tive chemical potential is then a good criteria to discuss 
the BCS-BEC crossover. 



IV. CONCLUSIONS 

A new type of density-dependent contact pairing inter- 
action was obtained to reproduce the pairing gaps in sym- 
metric and neutron matter obtained from a microscopic 
calculation [ll[. The contact interactions reproduce the 
two types of pairing gaps, i.e., the gap calculated with the 
bare interaction and the gap modified by medium polar- 
ization effects. It is shown that the medium polarization 
effects cannot be cast into the usual density power law 
form in symmetric nuclear matter so that another new 
isoscalar term g2 in Eq. ([9]) is then added to the density 
dependent term of the pairing interaction in Eq. ([T|) . 

We have applied these density-dependent pairing in- 
teractions to the study of the BCS-BEC crossover phe- 
nomenon in symmetric and asymmetric nuclear matter. 
We found that the spatial di-neutron correlation is strong 
in general in a wide range of low matter densities, up to 
kpn ^ 0.9 fm~^ [Pn/po ~ 0.15). This result is indepen- 
dent of the pairing interaction, either bare or screened 
one, as well as of the asymmetry of the uniform matter. 
Moreover, it is shown that the two pairing interactions 
mentioned above lead to different features for BCS-BEC 
phase transition in symmetric nuclear matter. To clar- 
ify the difference, we studied various order parameters, 
the correlation probability P((i„), the rms radius of the 
Cooper pair S^rms, the gap A„ and the effective chemical 
potential 

fm as a function of the Fermi momentum 
or equivalently as a function of the density. The screened 
interaction enhances the BCS-BEC crossover phenomena 
in symmetric matter, while the pairing correlations as 
well as the crossover phenomena are decreased in neu- 
tron matter by the medium polarization effects. For the 
screened-II interaction, the crossover reaches almost to 



the BEC phase at kpn ~ 0.2 fm"-'^ in symmetric matter. 
We should notice, however, that the BEC state is very 
sensible to the asymmetry of the medium and disappears 
in neutron matter. 
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APPENDIX A: EFFECTS OF THE CUTOFF 
PRESCRIPTION ON PAIRING GAP 

In the gap ([6]), the integral runs over the momentum fc, 
which is limited by the cutoff momenta k^ to avoid the 
ultraviolet divergence. There are several prescriptions for 
the cutoff momenta depending on the physical problem 
for which the interaction is applied. 

Prescription 1: This is the most simple prescription 
imposing on the single particle kinetic energy with 
the condition e„(fc) < Ec, i.e., A:+ — ^/2m*Ec/fi 
and k^ = 0. It is independent of the Fermi mo- 
mentum and of the pairing gap, but still has a 
weak dependence on the density through the ef- 
fective mass rn*ip) . It has been used by several au- 
thors [3, [H, and also adopted in shell model 
calculations in which all the shells up to a given 
cutoff energy are involved. 

Prescription 2: This prescription is based on the fact 
that the pairing occurs among states around the 
Fermi energy. Then, it is natural to define the 
cutoff energy with respect to the Fermi momen- 
tum by the condition e„(fc) < epn + Ec, i.e., 
k+ = ^2m*{€Fn + Ec)/h and stiU fc" = [1^. 
Through tp, the dependence on the density of this 
cutoff is much stronger than for the prescription 1. 
This prescription is close to the prescription used 
in HFB calculations. 

Prescription 3: This prescription is often used in the 
HFB calculations for which the cutoff is de- 
fined with respect to the quasi-particle energy 
^/{7Jk)^^lA^y^T^ < Ec. This leads to the fol- 
lowing definition of the cutoff momenta: k^ — 

1/2 



2m' 



A2 



/h (if Ec > A„). If 



k~ becomes imaginary for very small Vn, we set 
k~ = 0. In this prescription, the density depen- 
dence of the cutoff momenta kf^ is not trivial since 
it depends on the chemical potential and on the 
pairing gap. 
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E = 10 MeV E = 20 MeV E = 60 MeV 




kp„ [fm"'] kp^ [fm''] kp„ [fm''] 



FIG. 10: (Color online) The pairing gap of symmetric nuclear 
matter as a function of the neutron Fermi momentum kpn 
for the 3 prescriptions. The solid, dashed, and dotted lines 
are for the prescriptions 1, 2, and 3, respectively. Each panel 
corresponds to a fixed value for the cutoff energy, i?c=10, 20 
or 60 MeV. 



In Table |V| we give the parameters of the density- 
dependent term g = gi oi Eq. ([1]) to fit the bare gap using 
the prescriptions 1 and 2 for symmetric nuclear matter. 
Those for the prescription 3 are already given in Table [Til 
It shows how sensible these parameters are on the cutoff 
prescription. One could remark that the parameter 77^ 
is much more affected by the prescription than the pa- 
rameter as- We can also compare our results with other 
calculations. The parameters obtained with the prescrip- 
tion 1 can be compared to the one proposed in Ref. [la] , 
namely 77^ = 0.45, = 0.47 for Ec = 60 MeV. The value 
of the parameter as is significantly different. One can 
nevertheless obtain a comparable value for the parame- 
ter dig if one takes the approximation: Vn ~ ^Pn- With 
the prescription 2, we obtain similar parameters to those 
in Ref. [25] in which = 0.60 - 0.63, as = 0.55 - 0.58 
are obtained. The small differences can be explained by 
a different effective mass and different adopted pairing 
gap. 



APPENDIX B: COOPER PAIR WAVE 
FUNCTION 



It should be noted that in the limit Ec ^ A„, the 
prescription 2 and 3 give the same cutoff momenta . 
If the limit Ec ^ epn is also satisfied, then the three 
prescriptions are equivalent to each other. In Fig. fTOlare 
represented the pairing gaps obtained in symmetric nu- 
clear matter with the bare pairing interactions given in 
Table nil We compare the pairing gaps calculated for the 
three prescriptions with different values of the cutoff en- 
ergy, 10, 20 and 60 MeV. For a low value of the cutoff 
energy Ec=lQ MeV, the three prescriptions lead to very 
different pairing gaps. For the cutoff energy larger than 
20 MeV, the prescriptions 2 and 3 give very similar results 
because Ed ^1- Notice that at very low Fermi mo- 
mentum, kpn < 0.4 fm~^ {Pn/po < 0.01), the three pre- 
scriptions give a similar pairing gap, because both condi- 
tions Ec ^ A„ and Ec 3> cpn are satisfied. Nevertheless, 
for kpn > 0.4 fm~^, the prescription 1 gives a pairing 
gap different from the prescriptions 2 and 3, even for the 
larger cutoff energy Ec — 60 MeV. The reason is that the 
limit Ec 3> epn is not reached for kpn > 0.4 fm^"'^. 





Presc. 


1 


Presc. 2 


Ec [MeV] 




as 


ris as 


60 


0.461 


0.579 


0.593 0.537 


40 


0.413 


0.487 


0.657 0.506 



TABLE V: Parameters of the density-dependent term g ~ gi 
in Eq. ([l| obtained from the fit to the bare pairing gap in sym- 
metric nuclear matter. The effective mass is obtained from 
SLy4 Skyrme interaction. The parameters for the prescrip- 
tion 3 are already shown in Table HH 



From the Cooper pair wave function ()13|) in the mo- 
mentum space, we obtain the radial dependence of the 
Cooper pair wave function by a Fourier transform, 



pair 



(r) 



C 



(27r)3 
C I dk 



(fk UkVkC 

sin kr 

En{k) kr 



(Bl) 
(B2) 



where the normalization constant C is determined from 
the condition J drr^\^ pairir)]"^ = 1. The wave function 
'^pair is separated into two terms 4'i and '^2 defined in 
Eqs. (|16 p - (fT7)) . The term is solved numerically. If 
koo/ko 3> 1, with fco = \/2m* i^n/^i, we obtain the ana- 
lytical form (fT8|) for 'i'2{r) to the first order in M„/e„(fc), 
l/En{k) ^ 2m;jh^k'^. It is shown in Fig. [ID that the 
convergence is very fast and k^o = 2fco provides already 
a good converged solution. 

From the decomposition into ^'i and ^'2 presented in 
Eq. (fT5|l . the treatment of Matsuo is equivalent to set 
koo — kc and '^2 — 0, where kc is the cutoff momentum 
according to the prescription 2 in Appendix [XI By this 
treatment, a good agreement between the wave function 
obtained with the contact interaction and the one ob- 
tained with the Gogny interaction is reached for a cutoff 
energy Ec = 30 MeV. The agreement is very nice espe- 
cially in the region r < 3 fm. This result is easily un- 
derstood from the decomposition of 5* into 5*1 and 'I'2. 
Namely, by introducing a cutoff energy in the definition 
of the Cooper pair wave function, Matsuo's modification 
effectively removes the singularity at r ~ 0. The justifi- 
cation of cutoff in the Cooper pair wave function could be 
understood from model space considerations. For practi- 
cal reasons, the calculations for finite systems are never 
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FIG. 11: (Color online) The Cooper pair wave function r^'^paAr SiS well as its individual contribution r'^^'f and for different 
values of the ratio fcoo/fco and for two values of the neutron Fermi momentum. Notice that the convergence with respect to the 
ratio koo/ko is very fast. 



done with infinite basis but in a sub-basis involving a fi- 
nite number of wave functions. A cutoff is then naturally 
introduced in the finite model space and all the quanti- 
ties, including the Cooper pair wave function, are calcu- 
lated within the same model space. The cutoff treatment 
of Matsuo is thus common in nuclear matter and finite 
nucleus calculations. However, it should be remarked 



that despite the singularity of the wave function in nu- 
clear matter, the rms radius ^rms stays finite and larger 
than 3 fni, in the range of explored densities. We have 
checked that the treatment of Matsuo only affects the 
rms radius by about 10% in the worst situation, i.e., it 
increases the rms radius by less than 0.5 fm for neutron 
density at pn ^ Po/W. 
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